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Abstract 

The statistical thermodynamics of straight rigid rods of length k on triangular lattices 
was developed on a generalization in the spirit of the lattice-gas model and the classical 
Guggenheim-DiMarzio approximation. In this scheme, the Helmholtz free energy and its 
derivatives were written in terms of the order parameter <5, which characterizes the nematic 
phase occurring in the system at intermediate densities. Then, using the principle of minimum 
free energy with 8 as a parameter, the main adsorption properties were calculated. Compar- 
isons with Monte Carlo simulations and experimental data were performed in order to evaluate 
the reaches and limitations of the theoretical model. 



1. Introduction 

The adsorption of gases on solid surfaces is a topic of fundamental interest for various applica- 
tions. 1,2 From the theoretical point of view, the process can be described in terms of the lattice-gas 
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model.— — A lattice gas is a system of N molecules bound not more than one per site to a set 
of M equivalent, distinguishable, and independent sites, and without interactions between bound 
molecules. Many studies have been carried out on the adsorption behavior of small molecules in 
such systems. However, the problem in which a 2D lattice contains isolated points (vacancies) as 
well as £-mers (particles occupying k adjacent sites) has not been solved in closed form and still 
represents a major challenge in surface science. 

A previous paper 9 was devoted to the study of long straight rigid rods adsorbed on square lat- 
tices. In Ref. [— ], the Helmholtz free energy of the system and its derivatives were written in terms 
of the order parameter 8, which characterizes the nematic phase occurring in the system at inter- 
mediate densities i&i! Then, using the principle of minimum free energy with 5 as a parameter, 
the main adsorption properties were calculated. Comparisons with Monte Carlo (MC) simula- 
tions revealed that the new thermodynamic description was significantly better than the existing 
theoretical models developed to treat the polymer adsorption problem. 

In contrast to the statistic for the simple particles, where the arrangement of the adsorption sites 
in space is immaterial, the structure of lattice space plays such a fundamental role in determining 
the statistics of £-mers. Then, it is of interest and of value to inquire how a specific lattice structure 
influences the main thermodynamic properties of adsorbed polyatomics. In this sense, the aim of 
the present work is to extent the study in Ref. [-] to triangular lattices. The problem is not only of 
theoretical interest, but also has practical importance. A complete summary about adsorption on 
triangular lattices can be found in— - — and references therein. 

The rest of the paper is organized as follows. In Section 2, the theoretical formalism is pre- 
sented. Section 3 is devoted to describe the Monte Carlo simulation scheme. The analysis of the 
results and discussion are given in Section 4. Finally, the conclusions are drawn in Section 5. 
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2. Model and theory 



In this paper, the adsorption of straight rigid rods (or fc-mers) on triangular lattices is considered. 
The adsorbate molecules are assumed to be composed by k identical units in a linear array with 
constant bond length equal to the lattice constant a. The fc-mers can only adsorb flat on the surface 
occupying k lattice sites. The substrate is represented by a triangular lattice of M — L x L adsorption 
sites, with periodic boundary conditions. N particles are adsorbed on the substrate with 3 possible 
orientations along the principal axis of the array [see fc-mers marked with 1, 2 and 3 in Figure 1(a)]. 
The only interaction between different rods is hardcore exclusion: no site can be occupied by more 
than one k-mer unit. The surface coverage (or density) is defined as 6 — kN/M. 
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Figure 1: (a) Straight rigid rods adsorbed on triangular lattices. Solid circles (joined by thick 
lines) and empty circles represent tetramers (k = 4) and empty sites, respectively, (b) Schematic 
representation of the set of vectors N\ , N2 , N3 for a triangular lattice 



Let Ni, N2 and N3 be the number of rods oriented along directions 1, 2 and 3 on the surface, 
respectively. The total number of fc-mers is N = Ni + N2 + N3. According to DiMarzio's lattice 
theory (DiMarzio 1961), the number of ways Cl(No,Ni,N2,N3) to pack the molecules such that 
Af; of them lie in the direction i and there are A^o empty sites on the surface is given by 



a(M,N h N 2 ,N 3 ) = 



(i) 
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where N = M - kZj Ni. 

Since different fc-mers do not interact with each other, all configurations of N fc-mers on M 
sites are equally probable; henceforth, the canonical partition function Q(M,N\,N 2 ,N3, T) equals 
the total number of configurations, Q,(M,Ni,N 2 ,N3), times a Boltzmman factor including the total 
interaction energy between £-mers and lattice sites, ke^N 

Q(M,N u N 2 ,N 3 ,T)=q(T) N a(M,N u N 2 ,N 3 )ew(-Pk£ N), (2) 

where q(T) is the partition function for a single adsorbed molecule, /3 = l/kgT (being kg the 
Boltzmann constant and T the temperature) and e is the interaction energy between every unit 
forming a k-mex and the substrate. 

In the canonical ensemble the Helmholtz free energy F{M,N\,N 2l N 3 ,T) relates to Q. (M, N\ , N 2 , N 3 ) 
through 

PF(M,N h N 2 ,N3,T)=-lnQ(M,N h N 2 ,N3,T) 

(3) 

= -N\nq-\na(M,N u N 2 ,N 3 )+l3k£ N. 
Then, the remaining thermodynamic functions can be obtained from the general differential form— 

dF = -SdT - TldM + ^idN, (4) 

where S, U and /1 designate the entropy, spreading pressure and chemical potential respectively 
which by definition are, 

\dTj MN ' \dM J Jyy ^ \9 N J J jfi 
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Isotropic distribution of adsorbed &-mers 

For the case of an isotropic distribution of the £-mers, i.e., N\ = N2 = N3 = N/3, ?? reduces to, 



{[N + (2kN/3) + (N/3)]\y 
N Q \[(N/3)\} 3 (M\) 2 



(6) 



Applying the Stirling's approximation to ?? and replacing in ??, the Helmholtz free energy per site 
f = F/M can be written in terms of the intensive variables 6 and T, 



In 



+ %| + (i-e)in(i-e) 

k k 



(7) 



(0-3)ln3 + Se o 0--kWr). 

k 



Then, the chemical potential and the entropy per site s = S/M result 



jB/i =(*-!) In 



+ ln 



3£ 



k\n(\-6)-\nK e (T) 



(8) 



and 



5(0) 



In 



-ln--(l-0)ln(l-0) + (0-3)ln3 
k k 



e 



\nq(T)+T 



d\nq(T) 
dT 



(9) 



where K e (T) = q{T) exp (— /3&£ ) is the equilibrium constant. 

Anisotropic distribution of adsorbed A>mers 

To introduce the effect of the orientational order in the GD theory, it is convenient to rewrite the 
configurational factor in ?? in terms of the nematic order parameter <5,— 



(10) 
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8 represents a general order parameter measuring the orientation of the /c-mers on a lattice with m 
directions and the set of vectors {N\,N2, • • ■ ,N m } is characterized by the following properties: (z) 
each vector is associated to one of the m possible orientations (or directions) for a /c-mer on the 
lattice; (ii) the Nfs lie in a two-dimensional space (or are co-planar) and point radially outward 
from a given point P which is defined as coordinate origin; (Hi) the angle between two consecutive 
vectors, Ni and iVi+i, is equal to 2%/m; and (z'v) the magnitude of N t is equal to the number of 
/c-mers aligned along the /'-direction. Note that the iVi's have the same directions as the q vectors 
in.— These directions are not coincident with the allowed directions for the /c-mers on the real 
lattice. 

In the case of a triangular lattice, as studied here, m = 3, the angle between Nj and jvjf+i is 2n/3 
and ?? reduces to [see Figure 1(b)]: 

- N1+N2+N3 N1+N2+N3 ni , 

= = , (11) 

N1+N2+N3 N 

where \Nj\ = Ni has been used for notational convenience. 

5 can be expressed in Cartesian form as 5 = 8 x x+ Syf, where 

s = _J 2 2 £ (12) 

N V 

and 



In addition, 



kN = k(N { +N 2 +N 3 ) 
M M 
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Then, N\, N2 and N3 can be written as a function of 8 X , 8 y and 0, 



£4 (i-* + V5*). 

14 0-*-**). 



(15) 



Now, replacing ?? in the DiMarzio configurational factor ?? and using ??, the Helmholtz free 
energy per site can be written as, 



j8/(0,^,^)=(l-0)ln(l-0) 

(*-l) 



3fc 



(1 + 24)0 



In 



1-^-^(1 + 24)0 



1-^3^(1-^ + ^) e\ In [1-^^(1-4 + ^^)0 



1-^(1-4-^4) 



+ 



3k 

(1+24) 
3fc 







In 



1-^(1-4-^)0 



01n 



(1+24) 



3fc 







(1-4 + ^4) 

+ - 01n 



+ 



3fc 

(1-4-^4) 

3k 



9 In 



(1-4 + ^4) 

3fc 

(1-4-^4) 

3k 



+ j6e o 0--ln< ? (r). 



(16) 
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Finally, from ??, 



(k-1) 



and 



+ 
+ 



3 
3 

(*-l) 



1+2 



do 



o + 8 x ) 



In 



l-^r^(l+28 x )0 



3k 



1 -(V3 d -^ + d -^)9-8 x -V38 y 



V 30 90 



In 



3k 
3k 



1-5, 



1 

+ 3 



1+2 



+ 5,1 



In 



(1+25,) 
3fc 







+ 



+ 



(1-5, + ^) 




•ifcln(l-0)-ln^(r). 



h-V3 8y)j 







(17) 



s(e,5„5y) 
ks 



+ 
+ 



l-^^(l + 25,)0 
(*-l) 



In 



1- { \^(1 + 28 X )0 
3k 



(i-8 x +V3S^)e in i-t-11 (1-5^+^^)0 

In 1-^(1-5,-^)0 



01n 



(1+25,) 
3fc 

3fc 

(1-5 X -V35 y ) 
3fc 



(1 + 25,) 



-(l-0)ln(l-0) + 



3k 
01n 

01n 








(l-5, + y / 3 5 > .) 
3fc 

(1-5,-^5,) 



3k 



\nq(T) + T 











d\nq(T) 
~df 



(18) 



It is easy to see that, as |5| =0, i.e. 8 X = and 5y = 0, the isotropic case is recovered and, 
consequently, ?? reduce to ??-??. In general, the calculation of the adsorption isotherm and the 
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configurational entropy of the adlayer requires the knowledge of an analytical expression for the 
dependence of the nematic order parameter on the coverage. For this purpose, a free-energy- 
minimization approach can be applied. 9 The procedure is as follows: 

(1) We choose 8y = and 8 X ^ 0, this leaves N2 and N3 in an isotropic state, see ??. We can do 
this without losing any generality, since a pure nematic state is given by molecules aligned 
in one direction only. 18 

(2) By differentiating ?? (with 8y = 0) with respect to 8 X and setting the result equal to zero, the 
function 8(6) is obtained. 

(3) 8(6) is introduced in ?? and thus the adsorption isotherm and the configurational entropy of 
the adlayer are obtained (without orientational restrictions). 

The points (2) and (3) can be easily solved through a standard computing procedure; in our case, 
we used Maple software. 

3. Monte Carlo simulation 

In order to test the theory, an efficient hyper-parallel tempering Monte Carlo (HPTMC) simulation 
method— 1 ^ has been used. The HPTMC method consists in generating a compound system of R 
noninteracting replicas of the system under study. The i-th replica is associated with a chemical 
potential /!,-. To determine the set of chemical potentials, {fit}, the lowest chemical potential, 
jUi, is set in the isotropic phase where relaxation (correlation) time is expected to be very short 
and there exists only one minimum in the free energy space. On the other hand, the highest 
chemical potential, }1r, is set in the nematic phase whose properties we are interested in. Finally, 
the difference between two consecutive chemical potentials, \l[ and with /!,- < jU,+i, is set as 
A/j, = (jUi — jU/e) / (if — 1) (equally spaced chemical potentials). The parameters used in the present 
study were as follows: R = 25, fJL\ = — 10 and /j,r = 10. With these values of the chemical potential, 
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the corresponding values of the surface coverage varied from Q\ (jUi) ps 2 x 10 4 to 9r(i1r) 0.99 
for £ = 2, and from G^i) w 2 x 10~ 3 to 0*^) « 0.96 for it = 10. 

Under these conditions, the algorithm to carry out the simulation process is built on the basis 
of two major subroutines: replica-update and replica-exchange. 

Replica-update: The adsorption-desorption procedure is as follows: (1) One out of R replicas is 
randomly selected. (2) A linear fc-uple of nearest-neighbor sites, belonging to the replica selected 
in (1), is chosen at random. Then, if the k sites are empty, an attempt is made to deposit a rod with 
probability W = min [1 , exp(/3/i)] ; if the k sites are occupied by units belonging to the same /c-mer, 
an attempt is made to desorb this /c-mer with probability W = min[l,exp(— J3ju)]; and otherwise, 
the attempt is rejected. In addition, the displacement (diffusional relaxation) of adparticles to 
nearest-neighbor positions, by either jumps along the /c-mer axis or reptation by rotation around 
the /c-mer end, must be allowed in order to reach equilibrium in a reasonable time. 
Replica-exchange: Exchange of two configurations Xi and Xj, corresponding to the z'-th and j-th 
replicas, respectively, is tried and accepted with probability W = min [l,exp (—A)]. Where A in a 
nonthermal grand canonical ensemble is given by [— /3(/l/ — Hi) (Nj —Ni)}, and Ni (Nf) represents 
the number of particles of the z-th (j-th) replica. 

The complete simulation procedure is the following: (1) replica-update, (2) replica-exchange, 
and (3) repeat from step (1) RM times. This is the elementary step in the simulation process or 
Monte Carlo step (MCs). 

For each value of the chemical potential Hu the equilibrium state can be well reproduced after 
discarding the first ro MCs. Then, a set of r samples in thermal equilibrium is generated. The 
corresponding surface coverage 0/(jU,-) is obtained through simple averages over the r samples (r 
MCs). 

i (m) = -i,0[X i (t)}. (19) 

In the last equation, X, stands for the state of the z'-th replica (at chemical potential /i ; ). 

The configurational entropy S of the adsorbate cannot be directly computed. To calculate 
entropy, various methods have been developed. 21 Among them, the thermodynamic integration 
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method is one of the most widely used and practically applicable. The method in the grand canon- 
ical ensemble relies upon integration of the chemical potential ji on coverage along a reversible 
path between an arbitrary reference state and the desired state of the system. This calculation also 
requires the knowledge of the total energy U for each obtained coverage. Thus, for a system made 
of N particles on M lattice sites, 



In the present case U(N,M,T) = and the determination of the reference state, S(Nq,M,T), is 
trivial because S(Nq,M, T) = for A^o = 0. Then, using intensive variables, 



In this section, the main characteristics of the thermodynamic functions given in ?? will be ana- 
lyzed in comparison with simulation results and the main theoretical models developed to treat the 
k-mers adsorption problem. Three theories have been considered: the first is the well-known FH 
approximation for straight rigid rods;-^>21 the second is the GD approach for an isotropic distri- 
bution of admolecules;-2i>2£ and the third is the recently developed SE model for the adsorption of 
poly atomics 

The equations of the GD adsorption isotherm and the GD configurational entropy for an isotropic 
distribution of adsorbed rods were given in ??, respectively. The corresponding expressions in the 
FH and SE theories are as follows: 



S(N,M,T)=S(N ,M,T) + 




HdN' 



(20) 




(21) 



4. Results 





(22) 
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1- V , ; 



(25) 



The computational simulations have been developed for triangular Lx L lattices with L/k = 20 
and periodic boundary conditions. With this size of the lattice we verified that finite size effects are 
negligible. As mentioned in Ref. [ — ], the relaxation time increases very quickly as the £-mer size 
increases. Consequently, MC simulations for large adsorbates are very time consuming and may 
produce artifacts related to non-accurate equilibrium states. In order to discard this possibility, 
equilibration times ro of the order O(10 7 MCs) were used in this study. 

An extensive comparison among the new adsorption isotherm [??, solid line], the simulation 
data (symbols), and the isotherm equations obtained from the analytical approaches depicted as 
GD [??, dashed line], FH [??, dashed and dotted line], and SE [??, dotted line] is shown in Figure 
2: (a) k = 3, (b) k = 8 and (c) k = 10. In the case of??, 5(0) was obtained by following the 
minimization procedure described at the end of Sec. II. In addition, q(T) is set equal to one in the 
theoretical equations (vibrational and rotational degrees of freedom of the adsorbed molecules are 
not considered in the simulations). 

In part (a), the behavior of the different approaches can be explained as follows. The new 
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Figure 2: Adsorption isotherms for rigid fc-mers on a triangular lattice: (a) k = 3, (b) k = 8, and (c) 
k = 10. Symbols represent the MC results, and lines correspond to different theoretical approaches 
as indicated in part (a). The corresponding order parameters, obtained from the minimization of 
the free energy in ??, are shown in the insets, (d) Average percent error in the chemical potential, 
e^, as a function of k for the different approximations studied in this contribution 
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theory and GD agree very well with the simulation results for coverage values of up to 6 ~ 0.8; 
however, the disagreement between theoretical and simulation data increases for larger 6 values. 
The coincidence between the new theory and GD results is due to the fact that, for small values of 
k (k < 4), the function 8(6) minimizing the free energy is 8(6) =0 and, under this condition, ?? 
and ?? become identical. However, SE provides a good approximation with very small differences 
between simulation and theoretical results in all ranges of coverage. 

Let us consider now the case of k = 8 [Figure 2(b)] . The agreement between simulation and 
analytical data is very good for small values of coverage. However, as the surface coverage is 
increased, two different behaviors are observed. Although SE and the new theory provide good 
results, the classical FH and GD approximations fail to reproduce the simulation data. The differ- 
ences between GD and the theory in ?? are associated with the behavior of the order parameter 
8(6), which is shown in the inset of the figure. The functionality of 8 with coverage is indicative 
of the existence of nematic order for 6 > 0.4. Even though this result is not exact, the inclusion of 
8(6) in ?? leads to an extremely good approximation of the adsorption isotherm. 

The marked jump observed in the curve of the order parameter of a function of the coverage 
[see inset of Figure 2(b)] is indicative of the existence of a first-order phase transition in the adlayer. 
This behavior differs from that obtained for square lattices, 9 where the continuous variation of the 
order parameter with density indicates clearly the presence of a second-order phase transition in 
the adsorbed layer. This point is extensively discussed in the recent paper by. 28 

Figure 2(c) is devoted to the analysis of large adsorbates (k = 10, in the case of the figure). 
The results are very clear: (1) FH and GD predict a smaller 6 than the simulation data over the 
entire range of coverage; (2) SE agrees very well with the simulation results for small and high 
values of the coverage; however, the disagreement turns out to be significantly large in a wide range 
of coverage (0.3 < 6 < 0.9); and (3) in the case of the new isotherm, the results are excellent and 
represent a significant advance with respect to the existing development of k-mer thermodynamics. 

In order to compute the accuracy of each theory, the differences between theoretical and sim- 
ulation data can be very easily rationalized by using the average percent error in the chemical 
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potential e^, which is defined as, 

where /i s i m (/i appr ) represents the value of the chemical potential obtained by using the MC simu- 
lation (analytical approach). Each pair of values (jU s i m ,jU ap pr) is obtained at fixed 0. The sum runs 
over the N points of the simulation adsorption isotherm (in this case, N = 25 for all k). 

The dependence of on the k-mer size is shown in Figure 2(d) for the different theoretical 
approximations. Several conclusions can be drawn from the figure: 

1) In the FH and GD cases, increases monotonically with increasing k and the disagreement 
between MC and analytical data turns out to be very large (larger than 5%) for k > 5 and 
k > 6, respectively. 

2) For the SE theory, there exists a range of k (2 < k < 7) where £^ remains almost constant 
around 1.5% and SE provides a very good fitting of the simulation data. However, for k > 8, 
the differences between simulation and theoretical data increase with k. This deviation is 
associated with the appearance of an I-N phase transition in the adlayer for k > 7,— which 
is not covered by the SE theory. 

3) The agreement between the equation reported here [??] and the simulation data is excellent 
over the whole coverage range. This result provides valuable insight into how the adsorption 
process takes place. Namely, for k > 7 and intermediate densities, it is more favorable 
for the rods to align spontaneously because the resulting loss of orientational entropy is 
compensated for by the gain of translational entropy. 

4) The comparison with previous results obtained for square lattices 9 reveals that, for a fixed 
value of k, (i) increases with the connectivity for FH and GD theories, and (ii) in the case 
of SE approach and the new ??, does not change significantly as the lattice geometry is 
varied. 
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Mappr 



x 100%, 



(26) 



The differences between the approaches analyzed in this work can be also appreciated by com- 
paring the coverage dependence of the configurational entropy per site, which is presented in Fig- 
ure 3 for the same cases studied in Figure 2. and triangular lattices, respectively. The overall 
behavior of s(8) can be summarized as follows: for 0—^0 the entropy tends to zero. For low cov- 
erage, s(6) is an increasing function of 9, reaches a maximum at Q m , then decreases monotonically 
for 6 > 9 m . The position of m shifts to higher coverage as the fc-mer size is increased. In the limit 
0—7-1 the entropy tends to a finite value, which is associated with the different ways to arrange the 
fc-mers at full coverage. This value depends on k. 

As in Figure 2, GD and FH appear as good approximations in the low-surface coverage region, 
but the disagreement turns out to be significantly large for s(6 m ) and s(6 = 1). On the other hand, 
SE shows a good agreement with MC simulations up to adsorbate sizes of k ~ 8. Finally, in the 
case of ??, the agreement is notable for all 0, reproducing the MC results for s(6 m ) and s(6 = 1). 

As in the case of the chemical potential, an average percent error (e s ) was calculated for the 
difference between simulation and theoretical predictions. In this case, 



where s s i m (s a ppr) represents the value of the configurational entropy per site obtained by using the 
MC simulation (analytical approach). As in ??, each pair of values (Ssim^appr) is obtained at fixed 
6 and N = 25. 

The behavior of e s is similar to that observed in Figure 2(d). However, two main differences 
can be marked: (1) FH performs better than GD for all values of k, and (2) the differences between 
SE and ?? are more notorious, with ?? being the most accurate for all cases. 

Finally, analysis of experimental results have been carried out in order to test the applicability 
of the model proposed here. For this purpose, experimental adsorption isotherms of n-hexane in 
5A zeolites, previously compiled by Silva and Rodrigues, 29 were analyzed in terms of ??. Given 
that the experimental data were reported in adsorbed amount (g/100 g adsorbed) as a function 




(27) 
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Figure 3: Same as Figure 2 for the configurational entropy of the adlayer 
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of pressure, the theoretical isotherms were rewritten in terms of the pressure p and the adsorbed 
amount Q as fitting quantities. Thus, assuming that the adsorbed phase is in equilibrium with a 
ideal gas phase, the pressure p can be written asp« exp(/3/i). In addition, 6 = Q/Q max , where 
Qmax represents the maximum adsorbed amount. This choice allows us a direct comparison of ?? 
with the results obtained in Ref. [— ]. 
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Figure 4: Comparison between experimental and theoretical adsorption isotherms (adsorbed 
amount Q vs pressure p) for C^H^ adsorbed in 5A zeolite. Symbols represent experimental data 
from Ref. [ 32 ] and lines correspond to results from ??. The parameters used in the fitting procedure 
are listed in Table I 



As is common in the literature, 30 ' 31 a "bead segment" chain model of the molecules was 
adopted, in which each methyl (bead) group occupies one adsorption site on the surface. Un- 
der this consideration, k = 6 is set in the fitting data corresponding to C^. In this scheme, a set of 
isotherms of n-hexane in 5A zeolites for different temperatures were correlated by using only one 
value of Q m ax and a temperature dependent K e [T) as adjustable parameters. The results are pre- 
sented in Figure 4 and the fitting parameters are listed in Table I. A very good agreement between 
experimental and theoretical data is observed. In addition, the value obtained for the saturation 
adsorbed amount Q max = 12.1 is consistent with previous results reported in Refs. [— 

In summary, the analysis presented in Figures 2-4 demonstrates that (1) explicitly considering 
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the isotropic and nematic states occurring in the adlayer at different densities is crucial to under- 
standing the adsorption process of rigid rods, and (2) ?? provide a very good theoretical framework 
and compact equations to consistently interpret thermodynamic adsorption experiments of poly- 
atomic species. 

5. Conclusions 

The adsorption process of straight rigid rods of length k on triangular lattices has been studied 
via grand canonical Monte Carlo simulations, theory and analysis of experimental data. The pro- 
posed theoretical formalism, based on a generalization of the GD statistics, is capable of including 
the effects of the I-N phase transition occurring at intermediate densities on the thermodynamic 
functions of the system. 

The results obtained (1) represent a significant qualitative advance with respect to former de- 
velopments on fc-mer thermodynamics; (2) demonstrates that explicitly considering the isotropic 
and nematic states occurring in the adlayer at different densities is crucial to understanding the 
adsorption process of rigid rods; and (3) provide a very good theoretical framework and compact 
equations to consistently interpret thermodynamic adsorption experiments of polyatomic species. 
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Table 1 : Table of parameters used in the fitting of Figure 4 



Temperature (K) 


Qmax (g/100g flrf ,) 


K e (bar" 1 ) 


473 


12.1 


0.109 


523 


12.1 


0.402 


573 


12.1 


1.597 
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